SPECTRAL METHODS IN PDE 



W.-M. Wang 

Abstract. This is to review some recent progress in PDE. The emphasis is on (energy) 
supercritical nonlinear Schrodinger equations. The methods are applicable to other non- 
linear equations. 



1. Introduction 

We consider the nonlinear Schrodinger equation on the <i-torus T d = [0, 2-7r) d : 
d 

i—u = -Au+\u\ 2p u + H(x,u,u) (p>l,peN), (1.1) 
ot 

with periodic boundary conditions: u(t, x) = u(t, x + 2n7r), x G [0, 2ir) d for all n G Z d , 
where H(x,u,u) is analytic and has the expansion: 



H(x,u,u) = a m (x)\u\ 2p+2m u, 



with a m periodic and uniformly real analytic. The integer p in (1.1) is arbitrary. 

We study (1.1) from two perspectives, that of existence of quasi-periodic solutions 
and that of existence and uniqueness of smooth solutions to Cauchy problems. The 
link of these two questions are classical as the known invariant measure for smooth 
solutions are supported on KAM tori. Conceptually, one could view quasi-periodic 
solutions as providing a "basis" for the nonlinear solutions. 

It is well known [Bol] that on T d , when H = 0, (1.1) is locally wellposed in H s for 

s>max(0 (1.2) 
2 p 

which is derived by linearizing about the flow of the Laplacian and proving L p estimates 
of its eigenfunction solutions (Strichartz estimates). 

For d > 3 and sufficiently large p, the right side of (1.2) is larger than 1. For 
example, in dimension 4, the quintic nonlinear Schrodinger equation is locally well- 
posed in H s for s > 1, above the Hamiltonian H 1 topology, hence there is no available 
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conservation law. These equations are therefore supercritical as there is no a priori 
global existence from patching up local solutions, not even for small data, as (1.1) is 
non-dispersive, i. e., ||it||oo cannot tend to as t — > oo on the torus T d . 

The main purpose of this review is to describe the new method developed in [W2, 
3] to construct global or almost global solutions to (1.1) both in the perturbative and 
semi-classical regime. The semi-classical limit reveals further the geometric nature 
of this construction. This is not surprising, as already in the critical case, the local 
existence time depends on the Fourier geometry of the solution and not just its H 1 
norm. 

Following the PDE custom, when considering applications to the Cauchy problem, 
we will set the higher order terms H = 0. This is also because for small data, the 
construction carries over verbatim to H ^ 0. For KAM construction we choose to 
keep the translation invariance breaking term H, which traditionally lead to more 
complicated small divisors in the normalizing transform. Our method is different, it 
hinges on establishing a spectral gap, which is stable under small perturbations. So 
we keep H to illustrate this stability. 

To understand this stability, it is important to remark that the spectral gap is in 
the "space-time" sense as we put space and time on equal footing via a space-time 
Fourier series. 

The point of departure here is that this spectral gap is produced by the nonlinear- 
ity itself and not by eigenvalue variation of the linear operator, which has been the 
tradition. It is well known that eigenvalue variations are difficult to achieve with a 
multiplicative potential in dimensions 2 and above due to degeneracy of the Laplacian. 

A step was taken in [Wl] to break the degeneracy with a trigonometric polynomial 
potential in dimension 2, but much remains to be done. On the other hand, the spec- 
tral gap produced by the nonlinear ity is geometric in origin and hence robust. This 
is also because (1.1) essentially preserves translation invariance and hence the alge- 
braic properties of the exponentials, while adding a multiplicative potential destroys 
integrability. 

In the last part of the paper, we consider the nonlinear Schrodinger equation on gen- 
eral compact manifolds. We obtain critical Sobolev exponents for local well-posedness 
on general surfaces. This is an indication that the results on the flat torus may not be 
improved geometrically. 

2. Quasi-periodic solutions 

In this section, we consider the KAM aspect of (1.1). The Kolmogorov non- 
degeneracy conditions (or its weaker versions) are completely violated, as we perturb 
about a (infinite dimensional) linear system and not an integrable non-linear system. 
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Let be a solution to the linear equation: 



i"(o) = _ Au (o)_ (21) 
at 

We seek quasi-periodic solutions to (1.1) with b frequencies close to in the form 
of a nonlinear space-time Fourier series: 

u(t,x) = u{n,j)e in -° Jt e ij - x , (n,j) G Z b+d , (2.2) 

with u G M 6 to be determined. 

Writing in this form, a solution to (2.1) with b frequencies u/°) = {jl}k=i 
(jk 7^ 0) has Fourier support 

supp -u (0) = {(-e jk ,j k ),k = 1, ...,&}, 

where e Jfe is a unit vector in Z b and jk ^ jk' if k ^ k' . (Unless otherwise stated 
jk : = \jk\ 2 etc.) 

We note that (2.2) treats space and time on equal footing, contrary to the usual 
ODE approach, which views the solution as the time evolution of (Hamiltonian) vector 
fields. In this Hamiltonian language, the vector fields corresponding to the nonlinear 
Schrodinger equation in (1.1) are infinite dimensional. We also note that in (2.2), for 
each additional frequency in time, we add a dimension, which plays an important role 
in untangling the resonances. 

Define the bi-characteristics 

C = {(n,j) G Z b+d \ ± n ■ w (0) + j 2 = 0}. (2.3) 

C is the solution set in the form (2.2) to (2.1) and its complex conjugate in the Fourier 
space. This is the resonant set for the nonlinear equation (1.1). The fact that the 
singularities are not isolated points is the major difficulty here. 

To treat this manifold of singularities, we observe that in the absence of the non- 
linear term, there is no propagation in the Fourier space Z b+d : as there is no coupling 
among different (n,j), even though C is an infinite set. In the Fourier space, the non- 
linear term becomes a convolution operator and couples different (n,j). But as long 
as the connected (relative to this convolution) sets in C remain bounded in the first 
approximation, we should still be able to construct quasi-periodic solutions close to a 
given linear solution by implementing a Newton scheme. 

Let 

b 

u i0) (t,x) = J2a k e- ij ^e ijk - x , 

k=i 
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be a solution to the linear equation (2.1). We achieve this by first making a geometric 
selection in the Fourier frequencies {jk}k=i an d then an amplitude selection in the 
Fourier coefficients {cik}k=i- The geometric selection is in order to ensure bounded 
connected sets in C and is the main new ingredient in this theory. 

We call generic, if it satisfies the genericity conditions (i-iv) in [W2]. Here it 
suffices to say that the genericity conditions pertain entirely to the Fourier support of 
u ^'- {jk}k=i e and are determined by the |u| 2p w term in (1.1) only. Moreover 

the non-generic set O is of codimension 1 in (IR d ) b . 

The main result is 
Theorem 1 [W2]. Assume 

b 

fc=i 

a solution to the linear equation (2.1) is generic and a = {a k } G (0,5} b = B(0,S). 
There exist C, c > 0, such that for all e G (0,1), there exists do > and for all 
5 G (0, So) a Cantor set Q with 

meas {Q n B(0, 5)}/5 b > 1 - Ce c . 

For all a G Q , there is a quasi-periodic solution of b frequencies to the nonlinear 
Schrddinger equation (1.1): 

u(t, x) = J2 a k e- luJkt e lJk - x + 0(S 3 ), 

with basic frequencies u = {u>k} satisfying 

The remainder 0(S 3 ) is in an analytic norm about a strip of width 0(1) on f b+d . 

Remark. When d = p = 1, the non-generic set = 0. All are generic and 
only amplitude selection is necessary. This is the well understood scenario, see the 
Appendix in sect. 5. To understand the substance of the geometric and amplitude 
excisions in the theorem, it is useful to take H = and note the perpetual existence 
of periodic solutions: 

u = ae-^ + ^e^-* 

to (1.1) for all j G Z d and a G C. 

We have moreover the following semi-classical analog, which is new to the KAM 
context: 
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Corollary 1 [W2]. Set H = m fi.ij. Assume 

b 

k=i 

a solution to the linear equation (2.1) is generic, {jk\\ = i £ [KX d ] h , K G N + ane? 
a = {flit} £ (0, l] b = B(0, 1). There exist C , c > 0, such that for all e G (0, 1), there 
exists Kq > and for all K > Kq a Cantor set Q with 

meas {£n£(0,l)} > 1 - Ce c . 

For all a G Q , there is a quasi-periodic solution of b frequencies to the nonlinear 
Schrddinger equation (1.1): 

u(t,x) = J2 a ke~ l " kt e lJk ' x + 0(1/K 2 ), 
with basic frequencies u = {u>k} satisfying 

u k =j 2 + 0(l). 

The remainder 0(1/ K 2 ) is in an analytic norm about a strip of width 0(1) in t and 
0{1/K) in x on T b+d . 

Remark. These are quantitative, global, L 2 size 1 and large (kinetic) energy solutions, 
which could be relevant to the compressible Euler equations. 

A sketch of proof of the theorem 

We write (1.1) in the Fourier space, it becomes 

oo 

diag (n ■ u + j 2 )u + (u * v)* p * u + a m * (u * v)*^ p+m) *u = 0, (2.4) 

m=l 

where (n,j) G Z, b+d , v = u, u G M. b is to be determined and 

\a m (£)\<C'e- c 'W (C",c'>0) 

for all to. From now on we work with (2.4), for simplicity we drop the hat and write u 
for u and v for v etc. We seek solutions close to the linear solution of b frequencies, 
supp u(°> = {(—ej k ,jk),k = 1, ...,&}, with frequencies u/°) = {jk\k=i O'fc 0) and 
small amplitudes a = {ak} b k=1 satisfying ||a|| = 0(5) <C 1. 
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We complete (2.4) by writing the equation for the complex conjugate. So we have 



| diag (n-u+j 2 )u+(u*v)* p *u + J2™ =1 a rn *(u*v)< p+m ^*u = 0, 
I diag {—n ■ u + j 2 )v + (u * v )* p * v + J2m=i a ™ * ( u * i>)*^ +m ^ * v = 0, 

By supp, we will always mean the Fourier support, so we write supp for supp 
etc. Let 

S = supp U supp u^°\ 

Denote the left side of (2.5) by F(u,v). We make a Lyapunov- Schmidt decomposi- 
tion into the P-equations: 

F(u,v)\ Z b+ d \ S = 0, 

and the Q-equations: 

F{u,v)\ s = 0. 

We seek solutions such that u\s = u^°\ The P-equations are infinite dimensional and 
determine u in the complement of supp u^; the Q-equations are 2b dimensional and 
determine the frequency u = {^k}k=i- 

This Lyapunov-Schmidt method was introduced by Craig and Wayne [CW] to con- 
struct periodic solutions for the wave equation in one dimension. It was inspired by 
the multiscale analysis of Frohlich and Spencer [FS]. The construction was further de- 
veloped by Bourgain to embrace the full generality of quasi-periodic solutions and in 
arbitrary dimensions d [Bo4, 6]. More recently, Eliasson and Kuksin [EK] developed 
a KAM theory in the Schrodinger context. All the above results, however, pertain to 
parameter dependent tangentially non-resonant equations. 

We use a Newton scheme to solve the P-equations, with as the initial approx- 
imation. The major difference with [Bo4, 6, CW, EK] is that (2.5) is completely 
resonant and there are no parameters at this initial stage. The frequency a/ - 1 is an 
integer in Z & . So we need to proceed differently. 

First recall the formal scheme: the first correction 

A ($) = - CT = (2-6) 

where ^(i) ^ is the next approximation and F'(u^°\ v^) is the linearized operator 

on £ 2 (Z b+d x Z b+d ) 

F' = D + A, 

where 

D = ( dia § (n-^+j 2 ) 

y diag (— n ■ u + j 2 ] 
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and 



\p(u*v) p L *v*v (p + l)(u*v) p J ' v ' 

= A + O(S 2p+2 ). 

with u = uj(°\ u = -u(°) and v = v^°\ 

Since we look at small data, ||A|| = 0(5 2p ) <C 1 and the diagonal: ±n • u + j 2 
are integer valued, using the Schur complement reduction [SI, 2], the spectrum of F' 
around is equivalent to that of a reduced operator on £ 2 (C), where C is defined in 
(2.3) and to (D(5 2p+2 ) it is the same as the spectrum of A on £ 2 {C). 

To implement the Newton scheme using (2.5), we need to bound Aq 1 , which leads 
to genericity conditions (i-iv) in [W2]. From previous considerations, it suffices to 
consider Aq restricted to C. For generic u^°\ Aq\c = ®Ao, where Aq are Toplitz 
matrices of sizes at most (26 + d) x (26 + d). This can be seen as follows. 

Let 

C+ = {(n, j)\n- w<°> + j 2 = 0, 0} U {(n, 0) \n ■ = 0, m < 0}, 

C~ = {(n,j)\-n- u/ 0) + j 2 = 0, 0} U {(n, 0) |n • w<°> = 0, m > 0}, (2.7) 

c+n<r = 0, c+uc-=c. 



Assume (n, j) G C + is connected to (n',j') G C by the convolution operator Aq, then 
n' = n + An and j' = j + Aj, where (An, Aj) G supp (u^ * v^)* p , if (n', j') G C+ 
and 

(n-u/°)+j 2 )=0, 
(n + An) • + (j + Aj) 2 = 0; 

and if (n', j') G C _ , then (An, Aj) G supp («(°) * u' ')* 1 '- 1 * * «(°) and 

(n-u/°) +j 2 ) = 0, , 

(2 9) 

-(n + An) • u/°) + (j + Aj) 2 = 0. 

(Clearly the situation is similar if (n, j) G C~ .) 

(2.8, 2.9) define a system of polynomial equations. The generic condition on vS ' 
is imposed precisely to ensure that, aside from the degenerate case, (2.8, 2.9) have 
a solution for at most d + 2 equations. The degenerate case can be analyzed and 
correspond to systems of at most 26 equations reflecting translation invariance. So 
the largest connected set is of size at most max (26, d + 2) < 26 + d. Analysis of the 
degenerate case is in fact the reason for requiring the leading nonlinear 0(S 2p+1 ) term 
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in (1.1) to be independent of x. The x dependence of the higher order terms do not 
matter as they are treated as perturbations. 



The invertibility of A is then ensured by making an initial excision in a as is 
typically not an eigenvalue of a finite matrix. So x \\Aq 1 || < 0(5~ 2p ). Let 



F (u^,v^>) 

By requiring 



(0) fl (o h _/V°W >r*«< > 



supp F (« (0) y o) ) n{C\S} = 0, 
which is part of the genericity conditions in [W2], we obtain from (2.6) 

||Aw (1) || = ||At> (1) || < 0(S 3 ) 

for small 5. Inserting this into the Q-equations, which determine u, we achieve 
amplitude-frequency modulation: 

||Aa; (1) || x 0(5 2p ) 
|det(^— )| xO^" 1 ) >0 

ensuring transversality and moreover Diophantine u/ 1 ) on a set of a of positive measure. 
The tangentially non-resonant perturbation theory in [Bo4, 6] becomes available. 

Previously, quasi-periodic solutions were constructed using partial Birkhoff normal 
forms for the resonant Schrodinger equation in the presence of the cubic nonlinearity 
in dimensions 1 and 2 [GXY, KP]. These constructions rely on the specifics of the 
resonance geometry given by the cubic nonlinearity, see the Appendix in sect. 5. 



3. Almost global existence 

In this section, we set H = 0. As in sect. 2, we have results both on the perturbative 
as well as the semi-classical case. For the perturbative case, it is convenient to add a 
small parameter 5 and look at data of size 1. So we investigate the Cauchy problem 
on T d : 

f i-§^u = —Au + S\u\ 2p u (p > I,/? e N arbitrary), ^ 
\ u(t = 0) = Uq, 

Relying on the geometric information afforded by the resonance analysis in sect. 
2 and linearizing about a suitable approximate quasi-periodic solutions, we prove the 
following: 
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Theorem 2. Let u$ = u\ + 112- Assume u\ is generic satisfying (I. i-iv.) and \\u2W = 
O(S), where \\ ■ \\ is an analytic norm (about a strip of width 0(1) ) on T d . Let £(0, 1) = 
(0, l] b , where b is the dimension of the Fourier support of u\. Then there exists an 
open set A C £(0, 1) of positive measure, such that for all A > 1, there exists Sq > 0, 
such that for all 5 G (— 5q,5q), if G A, then (1.1) has a unique solution u(t) for 

\t\ < 5~ A satisfying u(t = 0) = uq and \\u(t)\\ < \\uq\\ + 0(5), moreover meas A — > 1 
as 5 ->■ 0. 

Remark. It is essential that the set A is open, which enables us to establish an open 
mapping theorem to analyze Cauchy problems . 

For perturbations of the Id cubic NLS (d = p = 1), similar stability results are 
proven in [Ba, Bo5]. For parameter dependent equations see [BG, Bo3]. The equations 
treated in [Ba, BG, Bo3, 5] are either L 2 or essentially L 2 well-posed. So there is a 
priori global existence. 

The equations treated in the theorem are of a different nature, there is no a priori 
global existence from conservation laws. In fact existence is obtained via explicit con- 
struction. Linearizing about approximate quasi-periodic solutions to prove existence 
and uniqueness for a time arbitrarily longer than local existence time, which is 0(d~ 1 ), 
is the main novelty. 

As in sect. 2, we also have the following semi-classical counterpart, providing quan- 
titative, almost global, L 2 size 1 and large (kinetic) energy solutions to Cauchy prob- 
lems. These solutions could be relevant to Cauchy problems for compressible Euler 
equations. 

Corollary 2. Set 5 = 1 in (3.1). Assume uq is generic with frequencies {jk}\ = i G 
[KZ d ] b , K G N+. Let £(0,1) = (0, l] b . Then there exists an open set A C £(0,1) 
of positive measure, such that for all A > 1, there exists Kq > 0, such that for all 
K > Kq, if {\uq\} G A, then (1.1) has a unique solution u(t) for \t\ < K A satisfying 
u(t = 0) = uq and ||w(t)|| < ||uo|| +C(l/-^ 2 )> where \\ ■ \\ is an analytic norm (about a 
strip of width 0(1/K) ) on T d , moreover meas A — > 1 as K — > 00. 

Remark. The previous related results are up to time 0{K), by solving the associated 
Hamilton- Jacobi equations before the arrival of caustics, cf. [Car] . 

A sketch of the proof 

Writing the first equation in (3.1) as F(u) = 0, for u satisfying the conditions in 
the theorem, we first find an approximate solution v such that 




(3.2) 
(3.3) 
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where r > A > 1. 

This approximate solution v is quasi-periodic with 0{\ log <5|) number of basic fre- 
quencies. The construction of v comprises of a finitely iterated Newton scheme to 
construct quasi-periodic solutions which solve (3.2) but (3.3) only to order 6 and then 
establishing an open mapping theorem using the spectral gap to solve (3.3) to order 
<V. 

Differentiate these quasi-periodic solutions with respect to the Fourier coefficients 
of it(°) gives a basis which spans L 2 (T d ) and allows to control the linearized flow. 
Schematically this could be understood as follows. 

Assume u is a solution satisfying the equation F(u) = and that it depends on a 
parameter a, then du/da is a solution to the linearized equation: 

f(«)( s ) = o. 

The difficulty here is that in the linearized equation for u, there is u. However the 
generic condition allows us to control the coupling of u and u on the flat torus. 

Using Duhamel's formula and the linearized flow to control the difference of (3.1) 
and (3.2, 3.3), we conclude the proof of the theorem. 

4. Nonlinear Schrodinger equations on surfaces 

We now consider the problem of nonlinear Schrodinger equations on smooth com- 
pact manifolds (without boundary): 

f i-§- t u = -Au + \u\ 2p u, 
\ u(t = 0) = uq. 

We specialize to dimension 2 and cubic nonlinearity (p = 1) as this is the most 
amenable. So we study the Cauchy problem 

f i§iu = —Au + \u\ 2 u, 
\ u(t = 0) = w , 

on a compact 2-manifold M. 

It is known from [BGT1, ST] that (4.1) is uniformly locally wellposed in H s for 
s > 1/2, in fact the flow is Lipshitz. We note that the corresponding Strichartz 
estimate on the flat 2-torus gives s > 0, cf. (1.2) and on the 2-sphere s > 1/4 [BGT2]. 
Both are known to be optimal [CCT, BGT2]. However the optimality of 1/2 on general 
surfaces M remained open. 

The following proposition constructs a counter example on the torus of revolution: 
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(4.1) 



Proposition. Let ds 2 = dx 2 + g(x)dy 2 with g G C 3 and admitting a unique global 
maximum. Then there are initial data, which are eigenf unctions of the Laplacian, such 
that the flow map is not Lipshitz in H s for s < 1/2. 

So the torus of revolution provides the geometric obstruction for going below H 1 / 2 
and 1/2 is the critical Sobolev exponents for local well-posedness on general surfaces. 

A sketch of the proof 

The Laplace-Beltrami operator decomposes into a direct sum of one dimensional 
Schrodinger operators indexed by the Fourier variable k in the y direction. We take 
as initial data the ground state ip of such a Schrodinger operator when k is large. Let 
A be its eigenvalue, the proof hinges on the equivalence estimate on the L°° norm of 
ip, namely 

IHUxA 1 / 8 (4.2) 

The proof of (4.2) uses Fourier and Hermite series analysis and a lemma of Bourgain 
[Bo2, D], which relates L 2 and L°° estimates. Afterwards one uses standard short time 
nonlinear analysis to complete the proof. (For the analysis details see [Cat].) 

We note that the above instability is in the sense of derivatives by exploiting the fact 
that in the linearized equation for u, there is u, and there is growth of L°° norm. On 
the flat torus, the exponentials have L p norm 1 for all p, which prevent this instablity. 

The above instability is uniquely infinite dimensional as it relies on non-equivalence 
of norms. However it is rooted in classical mechanics, namely the existence of lower di- 
mensional tori, which leads to eigenfunction concentration phenomenon. The Hermite 
series precisely represents the missing dimension in this picture. 

Remark. Using the same type of construction, 2/3 [BSS] should be optimal for surfaces 
with boundary. 

5. Appendix: the cubic nonlinearity on T d 

For simplicity we write u for and u for u/°\ the solutions and frequencies of 
the linear equation. The symbols of convolution for the cubic nonlinearity are \u\ 2 , 
u 2 and u 2 . Assume (n,j) G C + ((n,j) G C~ works similarly). In order that (n,j) is 
connected to (n'j') G C, it is necessary that either 

(a) [u*v](n,j;n'j') ^ or 

(b) [u*u](n,j;n'j') ^ 0. 
Case (a): Since 

n ■ u + j 2 = 0, 

n ' -uj + j' 2 = 0, 
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subtracting the two equations gives immediately 

Uk-j' k )-(j+jk) = 0, (5.1) 

where jk,jk' £ % d (k, k' = 1, b) and jk 7^ jk' if k ^ k', are the b Fourier components 
of u. 

Case (b): Since 

n ■ u + j 2 = 0, 
-ri ■ u + j' 2 = 0, 
adding the two equations gives immediately 

U + jk)-U + jk') = 0, (5.2) 

where jk,jk' £ ^ d (kj = 1; •••) 6) an d Jfc 7^ Jit' if ^ 7^ fc', are the b Fourier components 
of u. 

(5.1, 5.2) are precisely the well known resonant set for the partial Birkhoff normal 
form transform in [Bo4, GXY, KP]. (5.1, 5.2) describe rectangular type of geometry. 

supp F (u,v) n {C\S} = 

for the cubic nonlinearity in any d. When d = 1, (5.1, 5.2) reduce to a finite set of 2b 
lattice points in Z: {j = ±j k , k = 1, b} and O = in the theorems and corollaries. 
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